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Dartmouth College 

Semi-Annual Status Report - NCC 5-22 
February 1, 1980 - July 31, 1980 

This report covers activities of the Landsat Sensing Research 
Group (Earth Resources) and Of the Volcanic Gas Sensing Research 
Group (Planetary Science) which work in collaboration with the 
Goddard Institute for Space Studies, New York; Dr. Robert Jastrow, 
Director, Dr. Stephen Ungar of GISS is the Technical Officer for 
this project. 

NCC 5-22 supports work which was started in 1974 under NSG 5014. 

I. Landsat 

The Dartmouth Landsat Research Group continued application 
studies for Landsat data under the general category of analysis 
of vegetation cover, especially forestry and geobotany, that is, 
the effects of soil/earth mineral content on vegetation. 

• A. Forestry . The forestry group has had the same three 
permanent people working as in the previous year: Emily Bryant, 

Ken Sutherland, and Gibb Dodge. The'' were ably assisted by three 
undergraduates: John Stix, Dan Goodwin, and Dave Crampton. 

Our efforts lay in three general areas: Landsat classifica- 

tion research, software use and development, and spreading the 
word/contact with others. 

1. Classification Research 

a. John Stix and Em initiated a study of the use of the 
GISS fan algorithm in classification of forest types. The fan 
algorithm is used to classify a landscape consisting of a continuum 


between two pure cover types. In this case the continuum ranges 
from pure hardwood to pure softwood forest types. The fan algorithm 
places pixels in subcategories based on where their reflected radiance 
lies with respect to the signatures of the pure types. We wish to 
investigate the relationship between r+a fan subcategory that a 
pixel is placed in and the actual percent mixture of pure types 
found in the corresponding area on the ground. 

We have done background work on learning how to use the fan 
algorithm, getting the programs to run, and setting up the initial 
supervised classification from which the fan will select categories. 
When this is done, we will apply the classification to the Ashland 
District in northern Maine, an area which was classified previously 
in the Seven Islands project. 

b. Dave Crampton pursued the "fudge factor" question 
(temporal signature extension) using data from July and August of 
1976 covering the Ashland District in Maine. Using five sets of 
training sites, he made five estimates of a correction function 
(the "fudge factor"). It was applied to the July data as compensa- 
tion for atmospheric and other differences so that signatures used 
previously on the August data could be used on the July data as well. 

Preliminary results from classification of one township in the 
district show that the acreages of forest types in the July map 
differ from the August mao by 1.2% to 28% depending on the forest 
type and the correction function used. The function which was 
derived from all five sets of training sites combined gave acreage 
differences of -5.8% (softwood), +3.4% (mixed wood), and -4.7% 
(hardwood) . Further tests will be made using other estimates of 
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the function and applying it to another township. 

2. Software Development and Use. 

a. Dan Goodwin made trial runs of the various supervised 

and unsupervised classification algorithms in the ISURSL programs 
on two areas of interest to the forestry group: Concord, NH, and 

part of the Ashland District in Maine. The object was to become 
familiar with the programs, compare classification results with 
those from the GISS algorithm, and to identify the strong and weak 
points of the programs. 

He found use of the programs quite straightforward. Classifi- 
cation results compared well with the GISS algorithm— acreages for 
forest types on the Ashland district area were brought to within 
4.8% of a classification done using the GISS algorithm. Strong 
points in the ISURSL programs included a useful data display in the 
histogram program, the capability of running many training sites at 
a time to determine signatures, and the use of one of the programs 
for one-band density slicing (a kind of hand-made gray scale). 
Limitations encountered were that all jobs must be run from tape, 
there is no polygon masking program available, and that some of the 
statistics which are useful come with voluminous others which are 
not as useful and cannot be turned off. 

Me conclude from this project that classification results are 
dependent as much or more on the familiarity of the user with what- 
ever algorithm is being used as on the inherent attributes of the 
algorithm. 

b. Dan Goodwin also wrote programs to rotate classifica- 
tion map output, for instance, for the user who wants his digital 
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Landsat map to have north be vertical on the page. Ihese were 
written on both the Dartmouth and the GISS computers. The version 
on the GISS computer can handle a certain amount of scale change 
and change of symbols as well as rotation. 

c. Em wrote a program called DELTAS on the Dartmouth 
computer which helps the user choose classification parameters for 
the GISS algorithm. It displays training site data in "difference 
space” (rather than color space ) , making it easier to predict the 
effect that a change in classification parameters will have on 
classification results. The DELTAS program was applied to develop- 
ment of forest type signatures for the Ashland District. Comparison 
with accounts of an earlier effort at making signatures for the same 
area showed that the DELTAS program helped increase efficiency in 
signature development. 

d. Em wrote a program called ACCRACY on the Dartmouth 
computer. It displays a classification map with corresponding 
digital ground truth and prints a confusion matrix for the two. 

e. Signature packages which have been developed for 
forest types on various Landsat passes over the past five years 
were typed into the file EFESB . SIGPACK which is accessible through 
GISS/Wylbur . Any Wylbur user can use the signatures by copying 
them into a GISS classification program. 

3. Spreading the Word. 

a. The manuscript on the Seven Islands p iject was revised 
and was accepted by the journal Photogrammetic Engineering and Remote 
Sensing. (Abstract, Enclosure 1. A reprint will be included in a 
future report . ) 
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b. Landsat maps were given to John Ricard for field 
testing by local forest fire officials this summer. 

c. Help was given to students in Earth Sciences 32 
(remote sensing) in winter term with use of GISS hardward and 
software . 

1 d. Em gave a lecture on digital Landsat data classifica- 
tion to Dave Lindgren's (Geography) remote sensing class. 

\ 

e» Gibb and Em went to New York to talk over research 
plans with GISS personnel. 

f. Dartmouth remote sensing meetings have continued on 
the average of once a month. Sam Goward (Columbia) gave a guest 
talk in January. 

g. Maps of the White Mountain National Forest were 
presented to Forest Service personnel at their headquarters in 
Laconia, NH, by Gibb, Ken, and Em. Their preliminary reaction was 
that this type of information could be used better on a general 

level than on their local level, but we still await a more detailed 

* 

reaction. 

h. Gibb, Ken, and Em attended a meeting organized by the 

w 

Mew Hampshire State Office of Comprehensive Planning to investigate 
the possibility of setting up an ERRSAC Landsat Demonstration Project 
in New Hampshire. 

4 . Other items . 

We have had a fair number of hardware problems this half-year. 
The Cope card reader was down for a couple of week-long periods. 

Em will be entering the two-year Dartmouth Master's program 
in Computer and Information Science as a student in September and 


will be working half-time instead of full-time starting then. 

The progress report, furnished under the terms of the Dart- 
south Task Order to UNH, for the period October 31, 1979 - May 1, 

1980 is at Enclosure 2. 

B. Geology anu Geobotany . The geologic remote sensing group 
has worked in three major areas over the last six months. 

1. The project on the geologic mapping of the Landsat data 
has been completed. A Master's thesis was submitted by Joseph 
Francica (June 1980) (Enclosure 3) and a paper was presented at 
the 14th annual ERIM symposium in Costa Rica (Francica, Birnie, 
and Johnson, 1980) (Enclosure 4) . The project produced a series of 
It 100, 000 geologic maps interpreted from Landsat data. The maps 
cover over 11,600 square kilometers of the Ladakh Himalaya of 
Arabia and Pakistan. Landsat remote sensing techniques are nicely 
suited to this area because the hiqh relief and limited accessibil- 
ity make conventional geologic mapping very difficult and the dearth 
of vegetation and absence of atmospheric haze make clear bedrock 
exposures for the satellite. 

2. A project was undertaken using Landsat data to identify 

and map ‘lateral spectral contrasts in the Lower Hudson River Estuary . 
These lateral spectral contrasts were detected using the supervised 
classification algorithm developed at G1SS. The analysis of Landsat 
data lead us to conclude that there is a relatively strong down- 
stream flow on the west side of the river and that the western water 
is more turbid and less saline them the eastern water. These results 
were presented at the 14th annual ERIM Symposium in Costa Rica 
(Birnie and Posmentier, 1980) (Enclosure 5). 
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3. The third project involves our ongoing program to use 
airborne multispectral scanner data to detect geobotanical anomalies 
associated with mineralized ground. Data from the Mesatchee Creek 
experiment has been reevaluated and a paper has been revised and 
resubmitted to Economic Geology (Enclosure 6) . At Mesatchee 
Creek, a strong geobotanical anomaly was detected and was determined 
to correlate with the pyrite halo of a porphyry copper deposit. 

4. John Hughes wrote a number of programs for use with 
LANDSAT Digital Data on the Dartmouth Computer (Enclosure 7) . 

XI. Remote Sensing of Volcanic Emissions 
The principal activities of this group (Professor Richard 
E. Stoiber and Graduate Students Lawrence Malinconico and Stanley 
Williams) during this period (February - July 1980) were concerned 
with preparations, field work, and reports for visits to Guatemala 
and to Washington state (Mt. St. Helens). (Work in Nicaragua in 
February and August 1980 is related to this effort but is funded 
by the NSF. ) 

1. Guatemala. February 4 - 18, 1980. Stoiber, Malinconico, 
Williams. 

Our group was associated with the NCAR volcanic plume project 
in Guatemala. We flew through the gas of eruptive columns and 
beneath and through horizontal volcanic plumes. SO 2 output was 
measured with the airborne COSPEC. Excellent data were obtained 
from two volcanoes : Santiaguito and Feugo, (and some data from 

Pacaya) . It is being processed. An important feature of this 
project was that data were obtained by other scientists as to gas 


'composition and particulates. This was done at the same time as, 
and in coordination with, the SO 2 measurements. A discussion of 
the results is planned in Denver in late May. We will contribute. 

The new data collected chiefly concerned the configuration 
of the volcanic plumes , both in vertical and two horizontal 
dimension. Analyses such as these have not hitherto boen made. 

Also, see Enclosure 10, below, for discussion of review session. 

2. Mt. St. Helens. March 28 - April 6, 1980, April 9 - 13, 

1980, May 18 - 26, 1980, July 2-9, 1980. Stoiber, Williams, 
Malinconico. 

A major series of ventings and eruptions from Mt. St. Helens 
ever a period of several months (continuing) provided unique oppor- 
tunities to test and employ the COSPEC in an operational role. 

Four field trips, which were partially supported by NCC 5-22, 
are described in Enclosures 8, 9, 10, and 11. 

A fortunate confluence of interest with NBC news provided free use 
of an airplane for airborne monitoring and national publicity for 
the efforts and expertise of Professor Stoiber. 

In addition, the data obtained in the period March 28 - 
April 13, 1980, became the basis for an article submitted to SCIENCE 
on April 17, 1980 and published on June 13, 1980 (Enclosure 12). 

As reported in Enclosure 10, the Dartmouth-NASA COSPEC was 
loaned to the USGS on April 13, for use by them in continuous 
monitoring of Mt. St. Helens (until a USGS COSPEC could be made 
available.) Unfortunately, the scientist who was using the COSPEC 
was killed in the direct blast of the May 18 eruption and the COSPEC 
was destroyed, or at least buried irretrievably. Insurance coverage 


of the COSPEC (a common practice in field work using expensive 
instruments) provided funds for immediate replacement of the 
COSPEC, for use in subsequent work in Nicaragua. 

3. Nicaragua. February 1980, July - August 1980. Stoiber, 
Williams, Malinconico. 

This work, supported by the NSF, involves monitoring of the 
plumes of volcanoes Masaya and San Cristobal and generally helps 
develop our expertise in this field. 


Enclosures 1-12 a/s 


Landsat for Practical Forest Type Mapping 
A Test Case 

Emily Bryant 
Arthur G. Dodge, Jr. 

Samuel D. Warren 

Abstract : In a cooperative project, compute." classified 
Landsat maps were compared with a recent inventory of 
forest lands in northern Maine. Over the 196,000 
hectare (485,000 acre) area mapped, estimates of area 
of softwood, mixed wood, and hardwood forest types by 
the two methods agreed to within 5%. Cost of the 
Landsat maps is estimated at 6.5 cents per hectare (2.6 
cents per acre). Although the information derived from 
Landsat is not yet refined enough to he incorporated in 
current forest inventories, the techniques used are 
worth developing. 
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PROGRESS REPORT 



APPLYING LANDSAT MEASUREMENTS TO 
FOREST RESOURCE INVENTORIES 

October 31, 1979 - May 1, 1980 


Ken Sutherland and Gibb Dodge, Cooperative Extension Service, coordinated their 
activities with Emily Bryant, Dartmouth College, and other representatives of GISS. 


Site Selection 


1. Finished developing field maps of Groveton Papers (Diamond International) 
test compartment. Will deliver to the company and test usefulness. 

2. Assisted In developing criteria and categories for mapping White Mountain 
National Forest by districts. Met with WMNF Staff to obtain Input on maps developed. 

3. Developed map of the town of Plymouth for use by New Hampshire Forest Fire 
Service which was put Into the hands of local fire suppression personnel for test 
use. 


Guidance and Evaluation 

1. Guidance to: 

a) GISS on making changes In computer programs to produce better output 
products (polygon, legend, scale change, line, smoothing, etc.) con- 
tinued. 

b) Dartmouth work-study students - mapping techniques, observing ground 
truth sites, development of rotation and scale change program. 

-c) N.H. Forest Fire Service - using computer maps as field tools. 

. d) Revise and update 1980 work plans. 

2. Evaluation: 

a) Analyze computer outputs resulting from GISS program changes. 

b) Changes Indicated on Landsat maps by comparing 1975 and 1978 data. 

c) New rotation and scale changing program. 


Collaboration 


1. Serving on Cooperative Extension Service National Task Force for remote 
sensing - advice to Extension Committee on policy related to Extension activity 
In remote sensing technology transfer. 

2. Update information to the University of Vermont, University of New Hampshire, 
Milne Forestry Planning Group, GISS, ERRSAC, GSFC, Remote Sensing Group of Northern 


.New England, New Hampshire Fish and Gama Department and Office of State Planning. 

3. Remote sensing meetings with New Hampshire Office of State Planning, 
University of New Hampshire and Dartmouth. 

4. Net with N.H. Office of State Planning personnel on New Hampshire Landsat 
demonstration project with ERRSAC. 


Technology Transfer and Reporting 

1. landsat presentations to classes at Dartmouth. 

2. Generate remote sensing technology to New Hampshire Fish and Game Department 
and Division of Forests and Lands. 

3. Co-sponsor remote sensing workshop for State agencies with Office of State 
Planning. 


Submitted by. 


Kenneth I. Sutherland, Jr., 
Program Assistant 


Arthur G. Dodge, Jr., 
Program Leader 
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ABSTRACT 


Computer processed Landsat digital data were used to map 
the geology of 11,600 square kilometers of the Ladakh Himalaya 
of India and Pakistan. The results of a supervised 
classification algorithm were integrated with field studies to 
produce 1:100,000 geologic maps of this area. Six major rock 
types have been discriminated by Landsat data and were mapped. 
High relief and limited accessibility of the terrain make 
conventional geologic mapping difficult. This region was 
particularly suited for a Landsat data analysis, however, 
because the bedrock geology is clearly exposed as a result of 
the dearth of vegetation and the absence of atmospheric haze. 
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GEOLOGIC WI’HHfl Of THE LADAKH HIMALAYA 
BY COMPUTES PROCESSING OF LANDSAT DAT* 


Joseph M. Francica, Richard M. Birnie, and Gary D. Johnson 

Dept, of Earth Scioncca 
Dartmouth Colleqc 
Hanover, HH 0J75S 


ABSTRACT 

Coaq>uter processed Landsat digital data 
and field studies have been integrated 
to sake a geologic map of the Indus Su- 
ture in the Ladakh Himalaya. This co- 
ordinated approach has been successful 
at locating and identifying the areal 
extent of the major rock bodies in a 
2500 square kilometer area, much of 
which is inaccessable for conventional 
field geologic studies. 


1, INTRODUCTION 

Landsat digital data have bocn used to make a geologic map of portions of 
the Ladakh Himalaya of northern India and Pakistan. The study area encompasses 
the Indus suture, a tectonically complex tone of regional geologic importance. 
This region is particularly suited to Landsat-based analysis because of the 
clear exposure of the surface geology resulting from the paucity of vegetation 
and the absence of atmospheric haze. The high relief and limited accessibility 
of the terrain make geologic mapping by conventional techniques difficult. 

Approximately 12,000 square kilometers of terrain are being mapped; how- 
ever, only two representative regions totalling about 2500 square kilometers 
are presented here: 1) the Kargil region in the western limit of the study 

area, and 2) the Lch reqion in the eastern limit of the study area (Pig. 1) . 
These regions arc illustrative of our techniques and exhibit the major problems 
In the definition of rock body boundaries. 

2 . GEOLOGY 

The Indus suture, interpreted as a tectonic boundary between zones of 
varying rock type and petrogenesis (Heim and Gansser, 1939; Gansscr, 1964, 

1974, 1977, 1980; Powell and Conaghan, 1973 and 1975; Crawford, 1974) is char- 
acterized in the Kargil-Leh area as a zone of three to four narrow, steeply 
thrusted sheets of ophiolitic melange with associated ultramafic rocks (Frank 
•t al., 1977). The ophioli' .e melange and ultramafics, although well exposed 
and easily mapped in certain regions of the Kargil-Leh area, do not have 
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sufficient exposure to be mapped In the rcqions reported herein. The rock 
unite aeeociated with the Indue Suture and napped in thia study are described 

below . 


The Dras Volcanics (Wadia, 1937) are principally composed of metamorphosed 
eugita andesites, hornblende andesites and basalts (Trank et a i., 1977). In 
the eastorn portion of the study area, this unit grades into volcanoctastic and 
interbedded metasediments. Gupta et si., (197S) have established a Cretaceous 
age for this unit. 

The Ladakh Intrusivss (Frank ot jl,, 1977) have the largest areal extent 
of any rock unit within the study area. Generally, they vary in composition 
from granite to granodiorite in composition. Near Kargil, the intrusives in* 
elude a suite of norites, gabbros and diorites. To the cast near Lch, small 
diorite stocks and diabasic dikes are present. The intrusives are Cretaceous 
to Tertiary in age (Frank er al., 1977; Shah et al., 1976). 

The Ladakh Molassc (Towari, 1964) separates the Ladakh qranites to the 
north from the Dras Volcanics to the south. This late oroqcnic, primarily 
fluvial, clastic sequence is composed of intercalated channel sandstones and 
Ovarbenk mudstone lithofacies. 

The Indus Flysch (Frank ct al., 1977) consists of thinly bedded, fissile 
shales and sandstones. An early Nesoxoic age is suggested for these rocks 
(Towarl, 1964, Frank at al., 1977). 

Oaconaolidated Quaternary deposits abound in the Kargil-Leh area. 

3. LANDS AT DATA ANALYSIS 


The geocorrccted Landsat digital data (Scene E-301 35-04495) were classi- 
fied according to a supervised classification scheme developed by S. Ungar at 
NASA's Goddard Institute for Space Studies (Merry et al., 1977). This same 
classification scheme has been applied to a study of the Hudson River Estuary 
(Nirnie and Posmontier, this volume) where the following description is re- 
peated. In this classification scheme, a surface type is assigned a vector 
signature whose components arc the four individual Landsat band values. The 
spectral distance from this signature to a pixel to be classified is calculated 
and tested for acceptability relative to the user specified maximum allowable 
distance (DELMAX) . This dimensionless distance is based on two factors! color 
difference (AC) and brightness difference (AB) . 

To calculate the color difference (AC), the vectors in Landsat color apace 
corresponding to the signature and pixel to be classified are normalized to 
unit length. The AC is the chord length between the tips of the signature and 
piael unit vectors. 

The brightness of a pixel is determined by summing the four individual 
band reflected radiance values. The brightness difference (AB) between a pixel 
and signatura is defined as the absolute value of their difference in bright- 
ness. divided by their average brightness. 

The overall distance (D) between a pixel and a signature is calculated by 
the following relationship where the user chooses the degree to which color 
differences (AC) and brightness differences (AB) arc to be weighted (HCLR and 
NUT respectively) : 

D »>/mCLR» (AC) 1 ♦ WBKT’ (AB) ' 
where WBRT ♦ WCLR - 1. 


If the distance 0 it let* than a user-specified value (DELMAX) , then the un- 
' known pixel will be classified at belonging to the surface type associated with 
that signature. 

A four-dimensional classification volume is described about each signa- 
ture. The shape and orientation of the volume chanqcs depending on the values 
of WORT, MCLR, and DELMAX, but is generalised as a distorted hyper-ellipsoid 
inscribed within a truncated cone. 

Training sites of known gooloqic units were selected on tho basis of field 
work and published reconnaissance maps of the area. Spectral signatures for 
each unit were determined by averaqinq all pixels within a training site. Val- 
ues of WBRT, WCLR, and DELNAX for each signature were optimised by choosing 
those values that minimized the sum of the errors of omission and commission 
of each training site compared to all other training sites. 

The classification scheme is order dependant. Pixels which fall in more 
than one classification volume are classified at the first one tested. For 
example, areas of alluvium often classify as a granite because their volumes 
overlap and granite is classified before alluvium. A knowledge of which sur- 
face types created the overlap problems was used to set up the classification 
hierarchy and interpret the results. 

Reform field investigations were undertaken in August and September of 
1979, a preliminary map was generated for the Kargil region (Fig. 1). This 
area contains rock types representative of the Karqil-I.eh area as a whole. 

The town of Kargil served as a base of operations from which other parts of 
the area were accessible. Host of the training sites on which the signatures 
are based are located in this region. 

During the field investigations, training sites were studied in order to 
1) confirm their qcology and 2) identify the problems causinq incorrect 
classifications. T'.e accuracy of the classification of qeologic units outside 
the training sites was also studied at this time. These new observations led 
to the selection of new training sites, and a refined and improved version of 
the preliminary classification map was made (Figures 2 and 4). 

The large white areas represent unclassified pixels which, on the basis of 
geomorphology , we interpret to be shadowed. In these shadowed regions, the 
surface reflected radiance reaching Landsat is overwhelmed by backscattcrcd 
path radiance. The different reflectivity properties of various surface ma- 
terials are thus masked. Even though the Landsat data for Scene E- 101 35-04495 
(July 18, 1978) were chosen to minimize the shadow problems (sun elevation * 

58 degrees), the severe relief in the region results in extensive shadowinq. 

A geologic interpretation of the Landsat classification map was made by 
combining information in the classification map with field observations, qeo- 
norphic interpretations of 1:500,000 Landsat color composite imagery, and pre- 
viously published qeologic maps (refer to maps in Tewari, 1964; Shah rt al., 
197(i Frank at al., 1977). 

* 4 . DISCUSSION 


4.1 The Kargil Region 

The Dras Volcanics and the undifferentiated mafics of the Ladakh Intru- 
sive* have similar spectral signatures, respectively (Band 1 through Band 4 re- 
flected radiance values in sw , cm"’ , «r' li ch' 1 : Dras Volcanics * .607, .596, 

.SOS, 1.088; Mafics • .662, .596, .510, 1.1Z7). Because the volcanics appear 
first in our classification scheme, numerous mafic intrusive pixels are erro- 
neously classified as volcanics. This example of commission error is corrected 


by our field'lMicd obiorv.it ion th.it computer classified areas of volcanic;; ly- 
ing north of the main volcanic bolt and aeparated from it by other rock bodies 
■uat bo intarprotod aa mafic intruaivoa and not Draa Volcanlca. 

Tho Kargil baaln ha a boon the alto of fluvlo-qlaclal and lacuatrine sedl- 
awntation during tho lateat Tertiary and Ouatornary. Goomorphic evolution of 
tho Suru and Wakka river valleya during tho Quaternary haa reaultcd in tho de- 
velopment of a complex of terracea carved into theae baain-fill aedimenta and 
on certain of tho underlying older torrain. These wore rucognizcd on Landsat 
prior to entering the field but their areal extunt did not reconcile with the 
extant geologic mapa. Upon field study, we determined that the Landaat classl- 
f lea t ion accurately portrayed the diatribution of thia unit. Tho classifica- 
tion ahown in rigure 3 illuatratoa thia diatribution. 

Agricultural land, eaaily identified on Landaat color compoaite imagery 
and on our claaaif ication map, waa included aa alluvial depoeita becauae 
theoa lands are confined to irrigatible alluvial soils. 

The Boundary between the Draa Volcanlca and the Ladakh Nolaaso, which la 
clearly definod in the Landaat classification. aa a major rock body boundary, 
has been interpreted by Tewari (19641 as a thrust placing volcanlca over mo- 
laese. This thrust constitutes one of the imbricate faults associated with the 
Indus suture described above. 

4.2 The Lah Region 

The classification map for the Leh region shows clear outcrop patterns for 
both the molasse and the intrusives. Signatures used previously in the Kargil 
region for these units are applied here with equal success. A new signature 
was added for a molasse with a lower albedo. The molasse haa undergone a fa- 
cies change in this region and exhibits a different albedo when compared to the 
SMlaase of the Kargil region. 

A new signature for undifferentiated mafics of the Ladakh Intrusives un- 
successfully defined this rock type. The spectral signature for the mafics is 
similar to the signature for molasse (Band 1 through Band 4 reflected radiance 
values in mw*cm" ! 'sr* 1 "ch* 1 : Mafics - .824, .768, .634, 1.321; Molasse • .792, 

.729, .629, 1.264). The correction for the commission error in this region is 
accomplished by locating computer classified areas of molasse within the Ladakh 
Intrusives that we interpret to be mafics. 

The alluvial deposits in this region have the same spectral characteris- 
tics as the rock type from which they were derived. Alluvial fans south of 
the Indus River, derived from the darker molasso, classify as molasse. Tho 
same applies for fans derived from the Ladakh Intrusives north of the Indus. 
Their location and areal distribution are mapped by tho interpretation of hiqli 
oblique, ground-based photos whose orientations are known exactly. 

5. CONCLUSIONS 

The types of problems associated with the two test regions, Kargil and 
Leh, are common throughout the entire study area. Based on the criteria for 
Interpretation applied to these regions, other regions of the Ladakh Himalaya 
and adjacent Baltistan are being mapped in the same way. It is significant 
that -the signatures developed in the Kargil region were then successfully ap- 
plied in the Leh and other regions of the Indus suture zone. This fact illus- 
trates the application of computer processed Landsat digital data to geologic 
mapping of large areas. However, it should be emphasized that Landsat is no 
panacea. Only in conjunction with careful field observations can the problems 
associated with regional topography and geological complexities be resolved. 
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ARSTRACT 

An anatyaia of Landaat difital data haa allowed ua 
to da tact and map lateral apactral contraata in the Mud* 
aon River Catuary. Theaa eontraata ara believed to be 
related to turbidity in the iatuary with the waatarn aide 
taint more turbid. It ia clear that Landaat difital data 
ean be mapped and uaed aa an indicator of circulation and 
diaperaion in the Hudeon River Catuary. 


Is. INTRODUCTION 

A atudy haa been made of the Landaat difital data for the Lower Mud non Ri- 
ver Catuary (Scene C-109<-1$074, 21 Oct. 1*71*. The atudy area include a the 
ref ion of the Upper Say from i km aouth of the Verratano Narrowa up the river 
about. (S kb. to S km north of the Tappan See Sridfe (Pifure 1) . 

Initially, we are intereated in eonfirminf the preaence of lateral qra- 
dienta of watei proportion in the Hudaon River catuary to teat whether there 
are proceaaea that produce theae lateral fradienta at a rate auffieient to 
amiintain than in the preaence of lateral niainf. Until recently it haa been 
aaaiaaed that there are no lateral fradienta in the Hudaon River Catuary, but 
recent observational and theoretical data (Poementier and Raynent, ISIS) sup* 
port the preaence of lateral fradienta. 

'Secondly , if the lateral fradienta eaist, we are interested in determininf 
whether they produce a sufficient spectral contrast over a larqe enough area 
to make then detectable and nappable usinq Landaat difital data. 

2 . PRCVIOUS WORK 

Other inveatieators have uaed remotely sensed surface reflectance and re- 
flected radiance data to atudy eoastal and estuarine waters. Applications of 
remote senainq techniques to coastal tone nanaqement have been tummanted by 
Poaawntier (1*77). Rouse and Coleman (1S7S) correlated suspended sediment 
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MMMtttilM Md lilUctid NdiiMt in a laboratory experiawnt . They deter- 
mined that a spectral band equivalent to Undut MU load I <0.4-0. 7wnl gave 
tha but differentiation of suspended aedinent levels. Thoy assumed thia eor- 
ralatiaa oatandad to actual Land tat data and used band 1 Imagery to atudy tho 
alroulatloa of tha Louisiana bight. Johnaon and ioha (10??) found aiynlfleant 
oorrolatloaa between roflaotanoo in tho 0.440-0.4(0 and 0.430-0.440 m banda 
of aa atrborna null ispectral aoaanar and a nuobar of wator propartlr. measured 
dlrootly la tha Jama* River. Thaao properties included suspended aadiaont, 
chlorophyll a, nitrite, nitrate, phosphate, aalialty, and Secern depth. John- 
eon (1(71) found alallar reaulta for chlorophyll a and auaponded sediment in 
tho Mom York blyht. Khorran (1(7() uead four apoetral banda (0. 474-0. 4J0, 

D. (16-0. 410, 0.441-0.400, and 0. 744-0. 100 wal of aa Ocean Color dcanner on a 
0-1 aircraft to estimate aurfaee water electrical conductivity, ehlorophyll 
a c nc a a trotlon, auaponded eolide concentration, and turbidity. He we a able to 
■aa thoao four water propertiea in the ben rranciaeo bay Delta. Klemas and 
folia (1(74) uaed band S Landaat diyltal data and i angary to up auaponded and- 
iaoat c once ntration and tidal relatad varlatlona in auapondod aedinent distri- 
butiaa patterna in the Delaware bay. Landaat data were uaed to diatinyulah be- 
tween eeeweter and freeiMater, and aaoay four claaaea of a hallow marine and aa- 
teerlae water cateyorlea near Cape Pear (Nauael mt «i., 1(70). Maneaa and 
Mau a o l (1(74) delineated on the beaia of Landaat data diet they believed to bo 
tha aalt wedge boundary in the Pamlico Batuery . 

1. SALINITY DATA 

Aa part of a piiyalcal oceanography and interdiaclpllnary atudy of the 
Nmd oon batuery beginning in 1>71, temperature, aallnlty, and conductivity wore 
me aa u r ed at croaa-channel ttanaecta between the Apex of tha New York blg,:t to 
Oau g o cties. Mew York. Alorg each tranaoct, from two to four etatiena were oc- 
cupied! and the variablea at each a tat ion wera meaaured from the aurfaee to 
the bottom. Temperature and aallnlty meaturementa ware made with an Inter- 
O c a aa probe with a manually- logged digital readout and an analog plotter on 
dock, Tho seasonal, longitudinal, and vertical distrlbutiona of salinity and 
temperature have been analysed and are discussed by Posmentier and Raehlia 
(1(74) and Posmentier and Raynont (1(7(1 * 

Tho aallnlty data referred to in this paper were collected on 12-11 Aug- 
ust, H74, between 40*44* and 41*10' north latitude (Pigures 1 and 2). The 
times of observation span four semidiurnal tidal cycles. Any smali-acala sali- 
nity structure related to discontinuous point sources of seawater or fresh- 
water would, therefore, be beyond the resolving capacity of the survey. Polnt- 
by-point comparison of the non-synoptic salinities with synoptic reflected ra- 
diance data would be meaningless unless the latter were available for precisely 
the same tine as the salinity data, (inca time correlative Landaat data wera 
amt available, it was decided to com p are only the large-scale variations of 
aalialty with those of tha high guallty Landsat scene of 27 October, 1(72. 

taring the 12-11 Augusts 1(74 estuary study, there was very little temper- 
a turn variability in the aurfaee water. Therefore, only the salinity data, 
which waro measured with an accuracy of 0.04*/.,, are presented in Figure 2. 
There were fourteen individual latitudinal transects. The distance from tha 
shore to the ends of the transects was approximately 241 of tha Estuary width. 
Tho west aide surface salinities were lass than those on the east (Pleura 2) . 
This indicates that there is a greater proportion of freshwater on the west. 

The overage salinity difference was 0.740*/,.. An exception to the general 
tread occurred immediately south of the Piermont Pier (41*02' north latitude, 
Figure 1), where the higher salinity on the weat might be due to the large 
herlseatel eddy often observed on the wake of the pier which extends from the 
wasters shore one third of the way across "he Estuary. 
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There are aavaral factors which nay contribute to tha observed lateral 
salinity gradient in the Hudson Estuary. These include more freshwater runoff 
on the west side, tha Coriolis effect on the net downstream surface flow, the 
internal Kelvin waves forced by tides in the ocean, the curvature of the Estu- 
ary axis, asymmetry in tha Estuary cross section, and bottom friction. An 
evaluation of the hydrodynamic significance of these factors is beyond the 
scope cf this paper, but it is in preparation for publication elsewhere. What 
la eagjhasisad here is the interpretation that the higher salinity on the east 
aide of tha Estuary, whatever its cause, indicates - greater proportion of 
ocean water on the east side and a greater proportion of freshwater on the 
west side and confirms the presence of a lateral gradient of a water property. 

I 4. LANDS AT DATA ANALYSIS 

The geocorrected Landsat digital data (Scene E-1096-15074) were classified 
according to a supervised classification scheme developed by S. linger at NASA's 
Goddard Institute for Space Studies (Merry et •*!., 1977). This classification 
scheme has also been applied to a study of the distribution of rock types in 
the Ladakh Himalaya (Francica it <1. , this volume) where the following des- 
cription is repeated. In this classification scheme, a surface type is as- 
signed a vector signature whose components are the four individual Landsat 
•and values. The distance from this signature to a pixel to be classified is 
calculated and tested for acceptability relative to the user specified maximum 
allowable distance (DELMAX) . This dimensionless distance is based on two fac- 
torsi color difference (AC) and brightness difference (AB) . 

To calculate the color difference (AC) , the vectors in Landsat color space 
corresponding to the signature and pixel to He classified are normalized to 
unit length. The AC is the chord length between the tips of the signature 
and pixel unit vectors. 

Tha brightness of a pixel is determined by summing the four individual 
band reflected radiance values. The difference in brightness (AB) between a 
pixel and signature is defined as the absolute value of their difference in 
brightness divided by their average brightness. 

The overall distance (D) between a pixel and a signature is calculated by 
tha following relationship where the user chooses tha degree to which color 
differences (AC) and brightness differences (AB) are to be weighted (WCLR and 
NBRT respectively) s 

D -VWCLR«(AC) 1 ♦ WBRT* (AB) a 
where NBRT ♦ WCLR - 1. 

If the distance D is less than the user-specified value (DELMAX) , then the un- 
known pixel will be classified as belonging to the surface type associated 
with that signature. 

*A four-dimensional classification volume is described about each signa- 
ture. The shape and orientation of the volume changes depending on the values 
of NBRT, HCLR, and DELMAX, but is generalized as a distorted hyper-ellipsoid 
inscribed within a truncated cone. 

The initial step in the classification procedure was the selection of 22 
control areas in the Estuary. These areas measure 6x6 pixels and were cho- 
sen in east-west pairs at 11 latitudes (Figure 1) . T-tests were performed on 
the differences among the means of all 22 areas. A t-statistic of 1.692 or 
greater (95% confidence) was taken to indicate a significant difference between 
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two control trial . Two control ants Win choain on each aide of tha rivar (El 
and BE on tha aaat and Ml and W9 on tha west, Figure 1) auch that they wara si- 
■llar to most of tha othar control araaa on tha same alda of tha rivar and dif- 
farant fro* noat of tha control araaa on tha othar aide of tha river. Refer- 
ence signatures for waatarn and aaatern watar wara calculated by combining tha 
appropriate two control araaa for aach type. Landaat MSS Banda 4 ard 5 alona 
were sufficient to differentiate t'te waatarn and aaatarn refaranca aignaturaai 
therefore, Banda 6 and 7 wara diacarded, and a two-dinenaional classification 
wee performed . The raferanca aignaturea (aM.cm'' . ar" 1 .ch" 1 ) are liatad be low i 

li >i 

I western .446 .229 

aaatarn .294 .191 

Tha optimum values of NBRT, MCLR, and DELMAX wara chosen for aach refer- 
ence signature as tha lowest value of DELMAX and tha value of MBRT and MCLR 
that minimi tad tha sum of the errors of omission and commission. These values 
are liatad below and tha classification areas that they describe are shown in 
Figure 2 along with a B4-BS plot of tha pixels in tha two reference signature 
groups . 



DELMAX 

MBRT 

MCLR 

western 

O.OSS 

0.2 

0.9 

eastern 

0.05S 

O.S 

O.S 


Bach pixel in tha study area was subsequently claasifiad as *weatern*, 
■eastern*, or 'neither*. The classification is order dependent} once a pixel 
is successfully tasted for classification with one signature, it is not tested 
against other signatures. However, this was not a problem, because the over- 
all pattern and classification showed little difference whether ’eastern* or 
*weatern* watar was classified first. That is, there are very few pixels 
whose spectral characteristics placed them in both categories. Further, very 
few pixels in the Estuary fall in neither classification area and are unclassi- 
fied. Errors of omission and commission are given for the reference signature 
control areas in Figure 2. A copy of the printout of the pixel by pixel clas- 
sification of a representative portion of the river is shown in Figure 4. The 
boundary between ’eastern* and “western* water is interpreted as the point on 
each row where ’eastern* pixels begin to dominate over ’western* pixels. Un- 
classified pixels and classifications obviously influenced by ’sixth row strip- 
ing* are ignored. 


5. DISCUSSION OF THE LANDSAT CLASSIFICATION 


The interpreted classification map for the whole study area (Figure 1) 
confirms the existence of a lateral gradient in the surface spectral properties < 

of the Hudson River Estuary. The Landsat data clearly distinguish and map j 

these gradients on a scale at least as small as one half Kilometer. A t-test 
comparing the means of the reflected radiance values of the combined (11 x 26) 
eastern and combined (11 x 26) western control area pixels showed them to be 
significantly different for both Banda 4 and S with t-statistics of 12.44 and 
20.69 respectively. These are significant at greater than 99.95% confidence. 

Several possible explanations for this lateral gradient were mentioned earlier. ' 

Figure S shows a bathymetric cross section of the Estuary with the corres- 
ponding B4 and BS values of Landsat data plotted. Mhile the briqhter Landsat 
spectral valuas occur on the west where the river is shallow, it is clear that 
the Landsat spectral values do not increasa as the river shallows on the east- 
ern shore. This non-corralation of bathymetry and reflected radiance is also 

borne out on other cross sections of the Estuary not reproduced here. ) 
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The spectral properties of four freshwater reservoirs near the shore of 
the Hudson were also studied* and three of these freshwater reservoirs have 
spectral characteristics consistent with the more saline eastern river water. 
One of the reservoirs, Lake Tappan, had distinctly brighter, western water 
spectral characteristics. The water in the Lake Tappan reservoir was very 
silty because of ongoing coi.struction at the time of the Landsat overpass. It 
appears, therefore, that salinity is not the controlling factor in the spectral 
differences observed laterally across the Hudson, but that turbidity may con- 
trol the spectral differences. This observation is consistent with other stu- 
dies which were cited above and report increased reflected radiance values, 
particularly in Band S, with increased turbidity. 

There is no lateral contraat of Landsat spectral properties in the Estuary 
from 57th St. in Manhattan south to the lower portion of the Upper Bay. In 
this region, water with spectral properties of 'western* water predominated on 
both sides of the Estuary. This may be due to the .'act that the Landsat over- 
pass occurred at the time of slack following the ebb tide, which coincides with 
low tide. The tidal current chart (Figure 6) shows that this is the end of the 
(-hour period during which highly turbid East River and Kill Van Kull water and 
sewage effluent from lower Manhattan have been flowing into both sides of the 
Estuary. Specific locations of major sewage discharge into the Estuary are 
given by Simpson *e ml., (1975). 

Near the Verrazano Narrows, where very strong tidal currents and intense 
turbulence exist, the pixels are about evenly distributed between 'western* 
and 'eastern* classifications, with no discernable spatial clustering or pat- 
terns of either classification. 


6 . CONCLUSION 

This analysis of the Landsat data leads us to conclude that there is a 
relatively strong downstream flow on the west side of the river, and that this 
western water is more turbid and less saline than the eastern water. On the 
east side of the river, there is a weaker downstream flow, with saline water 
making its way farther upstream on the eastern shoreline. It is clear that 
differences in the spectral properties of the Landsat digital data can be 
mapped and used as indicators of circulation and dispersion in the Hudson Es- 
tuary. 
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Figure 1 . 

Lower Hudton Ri.er Estuary 


Interpretation of Landsat digital 
data for Lower Hudson River Estuary. 
Zones of generally "eastern" , "west- 
ern", and "mixed" water are ident- 
ified. Locations of control areas at 
eleven different latitudes are ind- 
icated as II through 111. At each 
latitude, a pair of control areas 
were chosen, one on the east and one 
on the west (for example, El and Ml 
respectively) . 
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rigure 2. Salinity Data 

Salinity as a (unction of latitude 
in tha Hudson Estuary. At each lat- 
itude, salinities for water on the 
eastern and western half of tha Est- 
uary are plotted individually. 
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Figure 3. B4-B5 Pixel Plot 

Plot of B4-B5 reflected radiance 
values for pixels in tha eastern 
(x) and western (+) reference sig- 
nature control areas. Values occur- 
ring in both control areas are 
shown with an s. The means of each 
control area are shown by the appro- 
priate bold lined symbol. Many 
values represent more than one pixel 
but these are not indicated. The 
classification area for each control 
area is outlined. 
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rigure 4. Classification 

Pixal by plxal classification of a 10.3 
kn stretch of tha Estuary between about 
40*56' and 41*02' north latitude. ■ rep- 
resents 'western* water, * represents 
'eastern* water, and blanks represent 
'neither* (unclassified) . The interpret- 
ed boundary between western and eastern 
water is drawn in. 
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ABSTRACT 


Visible and near infrared (450-1000 nm) reflected radiance 
spectra of the ground cover vegetation were measured at the 
Mesatchee Creek porphyry copper prospect in Washington. The 

p 

reflected radiance data are integrated over 400 m surfaces of 
a forest predominated by Douglas fir with lesser amounts of 
western larch. Analysis of the reflected radiance data indicates 
that the spectra from within the pyrite halo of the deposit have 
anomalously high reflected radiance values at 565 nm and low 
reflected radiance values at 465 nm. Six flight lines were flown 
on each of two days. Taking one flight line from the first day’s 
data as the control line, individual spectra with a 565 nm/465 nm 
reflected radiance ratio value greater than 1.7 fall predominantly 
within the p rite halo. When this threshold value is applied to 
all flight lines from the first day, 36.8# of the spectra within 
the pyrite halo and 5*4# of the spectra outside the pyrite halo are 
classified as anomalous. The zone of mineralization is clearly 

defined by the cluster of anomalous spectra. There is an 87# 

«* 

probability of an anomalous spectrum lying within the mineralized 
zone. The same technique was applied to the second day's data. 

The threshold ratio value was optimized at 1.6 and the resultant 
probability of an anomalous spectrum lying within the mineralized 
zone is then 92#. Higher cutoffs improve the probability of an 
anomalous spectrum falling in the mineralized zone; but fewer 
spectra are classified as anomalous, and the extent of the 



mineralized zone is not as well i t fined. The geobotanical anomaly 
correlates with the pyrite hale .nd is not preferentially con- 
centrated within the high Cu soil geochemical zones. Comparison 
with results obtained in lodgepole pine at Heddleston, Montana, 
show that different vegetation types manifest geobotanical anomal- 
ies in different spectral regions. 



JOHN M. HUGHES 
WINTER TERM, 1980 


The following computer programs were written for use 
with LANDSAT digital data: MATT, MATT1, TTEST, and BPLOT1. 
These programs are stored in R.W. Birnie's user number 
on the Dartmouth Computer. Printouts of the programs are 
attached. A brief description of each follows. 

MATT: Performs a student's T test on Control Area files 

(Pixel X Band files) . Outputs data in a file X file 
matrix, giving two places to right of decimal. Output 
represents confidence levels for the null hypothesis. 

MATT1: Same as above program, but prints out confidence levels 
in integer values. This enables user to print out 
up to a 24 X 24 matrix of control areas. 

TTEST: Program performs a t-test on only two different 

Pixel X Band files. Gives confidence levels for all 
four bands between two files. 

BPLOTL: Plotting program designed for use with Landsat 

MSS bands. The program asks for one or two files 
containing a Pixel X Band matrix and plots a 
point for each pixel in two band coordinates. 

Also calculates mean for each group of pixels, 
and draws Giss volume for each file if desired. 
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Trip Report 

Tot Mt. St. Helene, Vancouver, Washington 




On March 28, we flew to Portland, Oregon to witness 
the eruption of Mt. St. Helens, measure the SO^ output 
with the COSPEC, and obtain samples of fresh ash for further 
study. We secured the use of a Navajo Chieftain twin engine 
plane, courtesy of NBC News, for airborne use of the COSPEC. 

We had minor trouble with the sensitivity of the in- 
strument but were able to make 27 traverses under 11 different 
eruption plumes over the ensuring 7 days. SO^ output averaged 
approximately 30 tons per day. One fresh dry sample of ash 
was obtained. We returned to Hanover on April 6, 198Q. 

A manuscript, copy attached, was prepared, reviewing 
the scientific results of our visit and the subsequent visit 
by Larry Malinconico. It has been submitted to Science . 

The COSPEC, chart recorder, and tripod have been 
left in Vancouver, Washington. The COSPEC will be land 
carried to Denver before May 23. The accessories will be 
shipped to Dartmouth. David Johnston, of the USGS, now in 
Vancouver, WA. , is in charge of our instrument at this time. 


Richard E. Stoiber 
Stanley Williams 
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Dick Stoiber 
Stan Williams 
Phil Krueger 


Trip Report 

Tot Mt. St. Helens Volcano 
Vancover , Washington. 


On April 9, 1980 I travel 4, to Vancover, Washington 
with the COSPEC. For the next three days I worked with 
personnel from the U.S.G.S. measuring the S0 2 flux from 
Mt. St. Helens volcano. Several small eruptions were 
BMasured with flux rates averaging 30 t day -1 . On the 
11th and 12th a static plume was measured and had a 
flux of 3 t day -1 . 

I also instructed David Johnston and Tom Casadevall 
of the U.S.G.S. on the use and care of the COSPEC as they 
will be continuing the SO, monitoring of the volcano. 

On the 12th a sample of ash was collected from the 
summit of the volcano for leachate analysis. 

1 returned to Hanover on April 13, 1980. 


Lawrence L. Malinconico 
April 17, 1980 


EncJ 


TRIP REPORT 


MT. ST. HELENS AND NCAR REVIEW MEETING 

Stan Williams and Dick Stoiber spent May 18-22 and May 
25-26 at Vancouver, Washington. The purpose of the trip 
was to measure S0 2 after the big eruption of 18 May with 
our Cospec which we had left with Dave Johnston of the USGS. 
He. was making measurements without benefit of aiplane support. 
Since we knew NBC (Btfrrington) was at St. Helens we expected 
to be able to use our Cospec in a plane he would provide. 

On arrival in Vancouver WA. it was discovered that 
Dave Johnston had been killed by the direct blast of 18 May 
and our entire instrument package was destroyed. 

The USGS was expecting delivery of a new Cospec which was 
lost in transit. Upon arrival in Vancouver we turned our 
attention to tracing the instrument, expediting its delivery, 
arranging to send the Dartmouth chart recorder to Vancouver, 
and finding a portable 14 volt power source. We found the 
Cospec which had been broken in transit, and repaired it. 

We left for our NCAR meeting in Denver late Thursday 
by night flight from Portland, Ore. We went with Bill Rose, 
Irving Friedman, Leslie Zinser and Dick Cadle at NCAR. 
Publication and presentation of results of the Feb. 1980 
volcanic plume flights in Guat. were discussed with the 
following tentative plan: 

Abstracts for the Dec. 8-12 AGU in San Francisco as 


£nd IQ 


follows : 



1) Dispersion of SO 2 in volcanic plumes, Stoiber, 
Williams, and Malinconico 

2) SC >2 and isotopes at Pacaya and Santiaguito, 
Friedman, Williams, Malinconico 

1 3) SC >2 change through time at several volcanoes, 

Stoiber, Malinconico 

Publication of complete papers at near this date. 

We met with Barry Huebert of Colorado College on 
Saturday to discuss filters used to obtain HCl/SO^ ratio 
data. He will cooperate on gas work we intend to do in 
Nicaragua summer 1980 and written 1981. 

We left for Portland by night flight from Denver on 
Saturday, May 24. We monitored SO 2 from the air Sunday 
May 25. The SOj was being emitted from a medium sized 
eruption, which occurred near midday at a 2400 + day ^ 
rate. Data for further airborne measurements made while 
the volcano was relatively quiet the morning of May 26 has 
not been processed yet. 

We returned to Hanover Monday noon arriving early 
Tues. AMi( 

In summary: we find that SOj was now emitted at normal 
rates in the Mt. St. Helens eruption, we found the USGS 
instrument and put togethr an operative instrument package. 

We encouraged the USGS to use their instrument in the 

airborne mode and we think they will. We left our chart 

recorder with them til theirs is found (lost in transit) 

so measurements can be made. We collected an extensive suite 
of fresh ash for analysis in Hanover. We made a valuable contact 
with Barry Huebert with respect to future gas work. 
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